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(54) Fluoroscopy image reconstruction 



(57) A fluoroscopy imaging system (10) which, in 
one embodiment, includes a modified halfscan image 
reconstruction algorithm (34) which provides accepta- 
ble image quality along with the benefits of an enhanced 
temporal response as compared to the OS algorithm, is 
described. In an exemplary embodiment, the algorithm 



includes the steps of dividing the projection data into a 
preselected number of subsets, and for each subset, 
two partial images are generated: one with a first set of 
weights (denoted by and the other with a second, i. 
e., unity, set of weights (denoted by The final image, 
P h is generated from the two partial images. 
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Description 

[0001] This invention relates generally to imaging and more particularly, to a generating fluoroscopic images system. 
[0002] In at least one known imaging system generally referred to as a computed tomography (CT) system config- 

s uration, an x-ray source projects a fan-shaped beam which is collimated to lie within an X-Y plane of a Cartesian 
coordinate system and generally referred to as the "imaging plane". The x-ray beam passes through the object being 
imaged, such as a patient. The beam, after being attenuated by the object, impinges upon an array of radiation detec- 
tors. The intensity of the attenuated beam radiation received at the detector array is dependent upon the attenuation 
of the x-ray beam by the object Each detector element of the array produces a separate electrical signal that is a 

io measurement of the beam attenuation at the detector location. The attenuation measurements from all the detectors 
are acquired separately to produce a transmission profile. 

[0003] In known third generation CT systems, the x-ray source and the detector array are rotated with a gantry within 
the imaging plane and around the object to be imaged so that the angle at which the x-ray beam intersects the object 
constantly changes. A group of x-ray attenuation measurements, i.e., projection data, from the detector array at one 

is gantry angle is referred to as a "view". A "scan" of the object comprises a set of views made at different gantry angles, 
or view angles, during one revolution of the x-ray source and detector. In an axial scan, the projection data is processed 
to construct an image that corresponds to a two dimensional slice taken through the object. One method for recon- 
structing an image from a set of projection data is referred to in the art as the filtered back projection technique. This 
process converts the attenuation measurements from a scan into integers called "CT numbers" or "Hounsfield units", 

20 which are used to control the brightness of a corresponding pixel on a cathode ray tube display. 

[0004] To reduce the total scan time, a "helical" scan may be performed. To perform a "helical" scan, the patient is 
moved while the data for the prescribed number of slices is acquired. Such a system generates a single helix from a 
one fan beam helical scan. The helix mapped out by the fan beam yields projection data from which images in each 
prescribed slice may be reconstructed. 

25 [0005] In CT fluoroscopic systems ("CT Fluoro"), data collected from a helical or cine scan may be utilized to generate 
sequential frames of images to help, for example, in guiding a needle to a desired location within a patient. A frame 
corresponds to a two dimensional slice taken through the imaged object. Particularly projection data is processed at 
a frame rate to construct an image frame of the object. 

[0006] With known CT Fluoro systems, the general objective is to increase the frame rate while minimizing image 
30 degradation. Increasing the frame rate provides many advantages including, for example, that an operator physician 
is provided with increased information regarding the location of a biopsy needle. Typically, however, increasing the 
frame rate is at odds with minimizing image degradation. 

[0007] The time delay, or latency, of a CT Fluoro system is highly dependent on the nature of the reconstruction 
algorithm. For example, the time delay of one known CT Fluoro system is less with a halfscan image reconstruction 

35 algorithm (HS) and greater with an overscan image reconstruction algorithm (OS) . Even though the time delay is 
greater with OS algorithm, an OS algorithm may be selected in consideration of computational complexity. Specifically, 
since the weights utilized in the OS algorithm are channel independent, the order of filtering and weighting can be 
reversed. As a result, the data set can be divided into multiple subsets and two partial images can be generated with 
each subset. The final image can them be obtained by proper combination of various partial images. 

40 [0008] Such an approach, however, cannot be extended to the HS algorithm. Specifically, the halfscan weights are 
highly channel dependent. Therefore, the order of filtering and weighting can not be reversed. Also, the HS weighting 
function is also highly view angle dependent. Therefore, no advantage can be taken to make use of the previously 
filtered and backprojected data for generating a next image. If the halfscan algorithm can be implemented in similar 
efficiency to the above described OS algorithm, the temporal response of the CT Fluoro system can be significantly 

45 improved. 

[0009] These and other objects may be attained by a CT Fluoro system which, in one embodiment, includes a mod- 
ified halfscan image reconstruction algorithm which provides acceptable image quality along with the benefits of an 
enhanced temporal response as compared to the OS algorithm. In an exemplary embodiment, the algorithm includes 
the steps of dividing the projection data per gantry rotation into a preselected number of subsets, e.g., 8 subsets with 
50 each subset containing projection data which expands an angular range of 45 degree. Dividing the projection data into 
8 subsets enables, for example, achieving an 8 frames/s image update rate. For each subset, two partial images are 
generated: one with a first set of weights (denoted by 7)) and the other with a second, i.e., unity, set of weights (denoted 
by U t ). The final image, P,-, is generated in accordance with the following: 
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where N is the number of set-sets divided per gantry rotation from each data set. 

[001 0] An embodiment of the invention will now be described, by way of example, with reference to the accompanying 
io drawings, in which: 

[0011] Figure 1 is a pictorial view of a CT imaging system. 

[0012] Figure 2 is a block schematic diagram of the system illustrated in Figure 1 . 
[0013] Figure 3 is a radon space diagram for original halfscan. 
[0014] Figure 4 is a radon space diagram for a modified halfscan. 

15 [0015] Referring to Figures 1 and 2, a computed tomograph (CT) imaging system 10 is shown as including a gantry 
12 representative of a "third generation" CT scanner. Gantry 12 has an x-ray source 14 that projects a beam of x-rays 
16 toward a detector array 18 on the opposite side of gantry 12. Detector array 18 is formed by detector elements 20 
which together sense the projected x-rays that pass through a medical patient 22. Each detector element 20 produces 
an electrical signal that represents the intensity of an impinging x-ray beam and hence the attenuation of the beam as 

20 it passes through patient 22. During a scan to acquire x-ray projection data, gantry 12 and the components mounted 
thereon rotate about a center of rotation 24. 

[0016] Rotation of gantry 1 2 and the operation of x-ray source 14 are governed by a control mechanism 26 of CT 
system 1 0. Control mechanism 26 includes an x-ray controller 28 that provides power and timing signals to x-ray source 
14 and a gantry motor controller 30 that controls the rotational speed and position of gantry 12. A data acquisition 
25 system (DAS) 32 in control mechanism 26 samples analog data from detector elements 20 and converts the data to 
digital signals for subsequent processing. An image reconstructor 34 receives sampled and digitized x-ray data from 
DAS 32 and performs high speed image reconstruction. The reconstructed image is applied as an input to a computer 
36 which stores the image in a mass storage device 38. 

[0017] Computer 36 also receives commands and scanning parameters from an operator via console 40 that has a 
30 keyboard. An associated cathode ray tube display 42 allows the operator to observe the reconstructed image and other 

data from computer 36. The operator supplied commands and parameters are used by computer 36 to provide control 

signals and information to DAS 32, x-ray controller 28 and gantry motor controller 30. In addition, computer 36 operates 

a table motor controller 44 which controls a motorized table 46 to position patient 22 in gantry 12. Particularly, table 

46 moves portions of patient 22 through gantry opening 48. 
35 [0018] The image reconstruction algorithms described herein typically are implemented by image reconstructor 34. 

Such algorithms, however, could be implemented in other components of the imaging system such as in computer 36. 

Also, it should be understood that system 10 is described herein by way of example only, and the following described 

image reconstruction algorithms can be practiced in connection with many other types of imaging systems. 

[0019] Figure 3 depicts a Radon space representation of the HS sampling pattern. The lower and upper shaded 
40 triangles represent redundant sampling pairs. To ensure artifact free reconstruction, the redundant samplees must be 

properly weighted. An example of such a function is described below: 
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where y is the angle of a particular ray within the fan beam with respect to the central ray which intersects the iso- 
center of the CT system, and p is the gantry angle. 

[0020] Since the weights are highly channel dependent, the projection data has to be weighted by the weighting 
function before filtering. Also, the two shaded triangular regions occupy 180 degrees of projection view angular span. 
s That is, only about 20% of the entire projection data set can be switched for the order of operation for filtering and 
weighting. 

[0021] To overcome this difficulty, a modified halfscan algorithm (MHS) for the fast CT Fluoro reconstruction can be 
utilized. The Radon space for the MHS is segmented as shown in Figure 4. The two shaded triangular regions occupy 
only half of the original angular span of non-uniform weights. 
io [0022] The weights are: 



™ "V'/V jo ot h env i se 

(4) 

20 

[0023] At the lines p =y m -y, and p = % + y m -% the weighting function is not differentiate with respect to(. To suppress 
image artifacts caused by such discontinuity, feathering may be employed across the boundaries. A feathering distance 
of 20 channels is believed to be sufficient to suppress imaging artifact related to the discontinuity. 
The weighting function of the MHS has the following property: 

25 

tv(Y,p) = 1 - iv(y,tc + P) (5) 

Therefore, the weights for the upper triangle can be produced by subtracting the weights of the lower triangle from 
30 unity. That is, to produce a partial image with the weights of the upper triangle, two partial images can be produced: 
one weighted by the weights of the lower triangle and the other with unity weight. A difference image produced by 
these two partial images will be identical to the one produced directly from the upper triangle. 

[0024] In view of the above described weighting, the following image reconstruction algorithm can be utilized. Spe- 
cifically, the projection data set per gantry rotation is first divided into 8 subsets with each sub-set containing projection 
35 data which expands an angular range of 45 degree. Dividing the projection data into 8 subsets enables achieving an 
8 f rames/s image update rate. Alternatively, the data set could be divided into 60 degree subset to achieve 6 frames/ 
s image reconstruction rate, or other size subsets corresponding to a desired frame rate. 

[0025] For each subset, two partial images are generated: one with the lower (or upper) triangular weights (denoted 
by T 7 ) and the other with unity weights (denoted by Uj). The final image, P ; can be generated in accordance with: 

40 

(6) 

where N is the number of sebsets divided from each data set. 

[0026] The above described modified halfscan algorithm takes advantage of the symmetry in the weighting function 
50 and enables the reconstruction of images at 8 f rames/s with acceptable image quality. As a result, the modified halfscan 
algorithm enables the use of halfscan weighting in CT Fluoro applications to provide an image with acceptable image 
quality with the benefits of an enhanced temporal response as compared to the OS algorithm. 

55 Claims 

1. A fluoroscopy imaging system comprising a processor programmed to: 
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weight projection data in accordance with a ha if scan weighting function; and 

generate image data for partial images from the weighted data. 

s 2. A system in accordance with Claim 1 wherein to weight the projection data, said processor is programmed to divide 
the projection data set per gantry rotation into a predetermined number of subsets. 

3. A system in accordance with Claim 2 wherein the predetermined number of subsets corresponds to a desired 
frame rate. 

10 

4. A system in accordance with 2 wherein said predetermined number of subsets is equal to eight and each subset 
contains projection data which expands an angular range of 45 degrees. 

5. A system in accordance with Claim 2 wherein said predetermined number of subsets is equal to six and each 
is subset contains projection data which expands an angular range of 60 degrees. 

6. A system in accordance with Claim 1 wherein at least two partial images are generated, data for one partial image 
(Tj) being weighted with a first weight and data for another partial image (denoted by Uj) being weighted with a 
second weight. 

20 

7. A system in accordance with Claim 6 wherein said second weight is unity. 

8. A system in accordance with Claim 6 wherein a final image, P,- is generated in accordance with: 



25 



30 where N is the number of setsets divided from each data set per gantry rotation. 

9. A fluoroscopy imaging system in accordance with Claim 1 wherein said system is a computed tomography type 
imaging system. 

35 10. A computed tomography system for performing fluoroscopy scan, said system comprising an x-ray source, an x- 
ray detector aligned with said x-ray source, a patient table for being positioned between said x-ray source and 
said x-ray detector, and a display for displaying image reconstructed from data collected by said detector, said 
system configured to: 

40 weight projection data in accordance with a halfscan weighting function; and 

generate image data for partial images from the weighted data. 

11. A system in accordance with Claim 10 wherein to weight the projection data, said processor is programmed to 
45 divide the projection data set per gantry rotation into a predetermined number of subsets. 

12. A system in accordance with Claim 11 wherein the predetermined number of subsets corresponds to a desired 
frame rate. 



50 1 3. A system in accordance with Claim 1 0 wherein at least two partial images are generated, data for one partial image 
(T,) being weighted with a first weight and data for another partial image (denoted by U,) being weighted with a 
second weight. 

14. A system in accordance with Claim 13 wherein a final image, P,-is generated in accordance with: 

55 
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where N is the number of setsets divided from each data set. 

15. A method for generating an image using data collected in a fluoroscopy scan, said method comprising the steps of: 

weighting projection data in accordance with a halfscan weighting function; and 
generating image data for partial images from the weighted data. 

16. A method in accordance with Claim 15 wherein weighting the projection data comprises the step of dividing the 
projection data set per gantry rotation into a predetermined number of subsets. 

17. A method in accordance with Claim 16 wherein the predetermined number of subsets corresponds to a desired 
frame rate. 

18. A method in accordance with Claim 15 wherein at least two partial images are generated, data for one partial 
image (Tj) being weighted with a first weight and data for another partial image (denoted by U f ) being weighted 
with a second weight. 

19. A method in accordance with Claim 18 wherein a final image, P,-is generated in accordance with: 



where N is the number of setsets divided from each data set. 

20. A method in accordance with Claim 15 wherein the weighting function for weighting the projection satisfies the 
condition: 

= 1 - w(y,n + p) 
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